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The problem of finding the ground state of a frustration-free Hamiltonian carrying only two-body 
interactions between qubits is known to be solvable in polynomial time. It is also shown recently 
that, for any such Hamiltonian, there is always a ground state that is a product of single- or two-qubit 
states. However, it remains unclear whether the whole ground space is of any succinct structure. 
Here, we give a complete characterization of the ground space of any two-body frustration-free 
Hamiltonian of qubits. Namely, it is a span of tree tensor network states of the same tree structure. 
This characterization allows us to show that the problem of determining the ground state degeneracy 
is as hard as, but no harder than, its classical analog. 

PACS numbers: 03.67.Lx, 03.67.Mn, 75.10.Jm 



Quantum spin models are simplified physical models 
for real materials, but are believed to capture some of 
their key physical properties, which lie in the heart of 
modern condensed matter theory [1]. Ground states of 
strongly correlated spin systems is usually highly entan- 
gled, even if the system Hamiltonian carries only local 
interactions. So in general, finding the ground state of 
such a system is intractable with traditional techniques, 
such as mean field theory. 

In practical spin systems, different local terms in the 
Hamiltonian might also compete with each other, a phe- 
nomenon called frustration, which makes the system fur- 
ther difficult to analyze [2]. However, frustration is not 
a necessary factor to cause ground state entanglement. 
Frustration-free Hamiltonians can carry lots of interest- 
ing physics, ranging from gapped spin chains [3] to topo- 
logical orders [4, 5]. 

During recently years, the active frontier of quantum 
information science brings new tools to study quantum 
spin systems. In particular, local Hamiltonian prob- 
lems are shown to be in general very hard, i.e., QMA- 
complete [6]. It is also realized that the study of k- 
local frustration-free Hamiltonians for qubits is closely 
related to the quantum fc-satisfiability problem (Q-fc- 
SAT) [7], which is the quantum analogy of the classical 
^-satisfiability (fc-SAT), a problem that is of fundamental 
importance and has been extensively studied in theoret- 
ical computer science (see, e.g., [8]). 

Spin models with two-body interaction are of the most 
physical relevance, as two-body interaction, in particular 
of nearest neighbor or next nearest neighbors on certain 
type of lattices, are the strongest interaction terms in the 
real system Hamiltonian. Because two-level systems are 
most common in nature, spin-1/2 (qubit) systems are of 
particular importance. 



It is realized, however, that certain ground states of a 
two-body frustration-free (2BFF) Hamiltonian of qubits 
could be pretty trivial with almost no entanglement at 
all. Algorithmically, the problem of finding the ground 
state of a 2BFF Hamiltonian of qubits is known to be 
solvable in polynomial time [7]. It is also shown recently 
that for any such Hamiltonian, there is always a ground 
state that is a product of single- or two-qubit states; and 
if there is a genuine entangled ground state, the ground 
space must be degenerate [9]. There are also similar ob- 
servations of the ground states in random or generic in- 
stances [10-13], saying that the entire ground space is of 
a trivial structure, which is almost always the fully sym- 
metric space, with ground space degeneracy n + 1, where 
n is the number of qubits [10, 11, 14]. 

The main purpose of this work is to characterize the 
entire ground space in the most general setting. We im- 
prove the understanding of the ground space of 2BFF 
Hamiltonians of qubits by showing that it is always a 
span of tree tensor network states of the same tree struc- 
ture. In other words, these states are generated, from 
products of single qubit states, by the same series of 
isometries (from single qubit to two qubits). As this 
characterization holds for the most general case, it im- 
plies that computing the degeneracy of 2BFF Hamilto- 
nian (#Q-2-SAT) is in a complexity class called #P [15]. 
On the other hand, the classical analog #2-SAT of #Q-2- 
SAT is =^P-hard, therefore #P-complete.This answered 
a question raised in [11]. 

Two-body frustration-free Hamiltonian. — Consider a 
system of n qubits labeled by the set V = {1,2,. . . , n}. 
We will be interested in 2BFF Hamiltonians H Hj 
of the system. The Hamiltonian is called two-body if each 
term Hj acts non-trivially only on two qubits. The index 
J indicates the two qubits on which Hj acts. The Hamil- 
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tonian H is called frustration-free if its ground state also 
minimizes the energy of each term Hj simultaneously. 
Without loss of generality we can assume throughout 
the paper that the smallest eigenvalue of each term Hj 
is zero by shifting the energy spectrum. In this conven- 
tion, the frustration-free Hamiltonian H itself will have 
zero ground energy. Specifically we have 

1C{H) = f\(lC{Hj)®Hj), (1) 

where JC(H) is the ground space of H and T~Lj is the 
Hilbert space of the qubits not in J. From this equation, 
one easily sees that it is the ground space of each term 
Hj, not the structures of excited states, that matters for 
the ground space of a frustration-free Hamiltonian H. 
In other words, it suffices to consider local terms to be 
projections II j for our purpose. 

Closely related to the analysis of 2BFF qubit Hamil- 
tonians is the quantum 2-SAT problem (Q-2-SAT) first 
considered by Bravyi [7] . Naturally generalizing classical 
2-SAT, the Q-2-SAT problem asks whether, for a given 
set of two-qubit projections {IIj} of an n-qubit system, 
there is a global state |^) such that iljl^) = for all J. 
Apparently, we answer "yes" to the problem if and only if 
the Hamiltonian ^2 Hj is frustration-free. It was known 
that Q-2-SAT is decidable in polynomial time on a clas- 
sical computer [7] . The proof of the statement actually 
constructs a specific n-qubit state in the ground space 
of n j if there is any. Our techniques will be similar to 
those used by Bravyi, but we will show a stronger result 
that one can not only find one state in the ground space, 
but also represent the entire ground space in terms of a 
span of special states. 

A case study of the rank. — Given a 2BFF Hamiltonian 
H = Y^Hj, what can we say about the ground space 
JC(H)1 First of all, as argued previously, we only need to 
consider Hamiltonians of the form H = ^ II j where nj's 
are projections onto JC{Hj)^~. We will start our analysis 
by considering the rank of the projections II j. 

First, if there is a Uj of rank 3, the only possible state 
for the two qubits in J is / — Hj of rank 1, and this 
reduces to a problem on qubits in V \ J. 

If there is a Hj of rank 2, the state of qubits in J 
is restricted to a two-dimensional subspace. Let 1^0)0,6 
and \ipi) a ,b be two orthogonal states that span the sub- 
space, where a, b are the two qubits in J. One can en- 
code qubits a and b by a single qubit d. For this pur- 
pose, we define an isometry U in the following form 
U : \0)d >-> \i>o) a ,b, |l)d ^ \i>i)a,b- This procedure pro- 
duces a set of constraints on n — 1 qubits. It is easy to 
verify that a state is in the ground space of the re- 
duced problem if and only if U^l!) is in the ground space 
of the original problem [7, 9]. 

When there is no projection of rank larger than 1, we 
are dealing with the homogeneous case [7]. It turns out 
that the homogeneous case is the hardest and we will dis- 
cuss it two separate sections. As we will see, the ground 



space of the homogeneous Hamiltonian (more precisely, 
the simplified homogeneous Hamiltonian defined later) is 
spanned by single-qubit product states. The above case 
analysis gives an explicit representation of the ground 
space of a general 2BFF qubit Hamiltonian, which is 
given by the following 

Main Observation — The ground space is always a 
span of tree tensor network states of the same tree struc- 
ture. 

We illustrate this observation in Fig. 1, where the 
ground space is viewed as a span of states generated 
by the isometries (blue triangles) organized in a forest 
form (a collection of trees) acting on product states (in- 
put from the left). In the language of tensor network 
states [16, 17], one can also represent these states in terms 
of tree tensor network after combining the input product 
states and the roots of trees in the forest. 



Input: a set of product 
states that spans the ground 
space of the simplified ho- 
mogenous Hamiltonian. 




FIG. 1: The general structure of the ground space 

Homogeneous case with product constraints. — Con- 
sider the Hamiltonian H = J2 Tlj where Hj's are rank-1 
projections. One can visualize the interactions in H by 
a graph G. The graph has n vertices corresponding to 
the qubits and two vertices are connected when there is 
a non-trivial interaction Hj acting on them. We will also 
distinguish two types of edges in the interaction graph. 
Let n = \<fi)((f>\ be a projection. We will use a solid edge 
in the graph when |0) is entangled and a dashed edge 
when \4>) is a product state. Let us first focus on the 
homogeneous case with product constraints only. 

In this case, the interaction graph consists of dashed 
edges. We will show that the ground space is a span 
of product of single-qubit states (or, for simplicity, a 
product span). It will also be useful to know that the 
states we choose are orthogonal up to a local operation 
L = <2)j=i Lj, where Lj is a non-singular local operator 
acting on the j-th qubit. Note applying L on the 2BFF 
Hamiltonian H = J^Hj results in H L = Y^Lj HjLj, 
where Lj = QQ^jLj. And H L , which is also 2BFF, 
has the same ground state degeneracy as H [7, !)]. The 
relation between the ground space of H and H L is 

L^JCiH) =JC{H L ). (2) 

Before we actually give the proof, let us first examine 
several simple examples. The first example considers a 
chain of interactions as in Fig. 2b. Let \ctj) <g> \ j3j) be the 
constraint on the j-th edge. We will call it an alternating 
chain if |/3j-i) and \oij) are linearly independent for all 
j. It is easy to see that the solution space is k + 1 for 
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an alternating chain of k qubits. The second example 
shown in Fig. 2c is called the alternating loop. As its 
name suggests, it is a loop where the two constraints 
on any vertex are linearly independent. Any alternating 
loop has solution space of dimension 2, namely the span 
of 1 00 . . . 0) and |11 . . . 1) up to the local operation that 
maps | ay) and to |0) and |1). The final example we 

consider is called the quasi-alternating loop. It is almost 
the same as the alternating loop except that there is one 
special vertex on the loop having the same constraint on 
the two edges adjacent to it. Figure 2d gives such an 
example where the top vertex is special. It is easy to 
see that the constraint on the special vertex of a quasi- 
alternating loop must be satisfied. In particular, for the 
loop in Fig. 2d, the top vertex must be |1) as otherwise 
it will be impossible to satisfy all five constraints on the 
loop. 
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FIG. 2: Dashed interaction graph and three examples 

We now start the proof by induction on n, the number 
of qubits. For n = 1,2, the observation is trivial. If 
there is a vertex a on which the constraints are the same 
up to global phases, let the constraints be |0) a and, more 
concretely, let the constraints on an edge connects to a be 
of the form |0) a |a)b for some qubit b. We can write any 
state in the ground space as |$) = |0) a |$o) + |l) a |3>i)- 
Obviously, |$o) an d l^i) are both in the ground space of 
the constraints not acting on a. Moreover, |$n) also needs 
to be orthogonal to |a)(,'s. By the induction hypothesis, 
both |$ ) and |<3?i) are in a product span. Therefore, 
|$) is also in a product span. On the other hand, if one 
cannot find any vertex whose constraints are the same, we 
can find either an alternating loop or a quasi-alternating 
loop in the graph. If a quasi-alternating loop is found, we 
know the state for the special vertex of the loop and can 
use the induction hypothesis on the remaining system. 
Otherwise, if an alternating loop is found, we can write 
any state in the ground space as 



|$) = |00 -..0)1*0} 



111 - - - 1> l^i 



(3) 



up to local operations on the loop. If a constraint acts on 
two qubits on the loop, it can only restricts the loop to 
be exactly |00 • ■ ■ 0) or 1 1 1 • ■ ■ 1). The analysis is similar 
to the first case when a constraint I a) a 1/3) {, acts on one 



qubit a on the loop and another qubit b outside of the 
loop. This completes the proof. Notice that the local op- 
erations chosen here are determined by the constraints of 
alternating loops, and that one will never have two alter- 
nating loops giving different local operations for a single 
qubit, the orthogonality of the states up to local opera- 
tions follows. We note that the orthogonality property 
only holds for the product constraints. The symmetric 
subspace, for example, is not a span of orthogonal prod- 
uct states up to local operations although it is the span 
of |00), |11), |++) where |+) = (|0) + \l))/y/2. 

General homogeneous case. — Given a general homo- 
geneous Hamiltonian, the interaction graph will consist 
both solid and dashed edges. The main technique is to 
simplify the interaction graph in hand without changing 
the ground space. Two sliding operations as shown in 
Figs. 3a and 3b will be used in the simplification. The 
Type -I sliding says that if we have entangled interactions 
between 1,2 and 1,3, we can change it to two entan- 
gled interactions between 1,2 and 2,3 without affecting 
the ground space. The Type-II sliding is of a similar 
spirit, but involves both entangled and product interac- 
tions. We will only prove the validity of Type -I sliding 
as a similar argument holds for the Type-II sliding op- 
eration. 
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FIG. 3: Simplification of the interaction graph 

Let II12 = |0) (0| and II13 = IV'XV'I be the two rank-1 
operators acting on qubit 1,2 and 1,3. We will find a 
local interaction II23 acting on 2,3 such that II12 + II23 
has the same ground space as II12 + II13. As \<j>) and 
\ip) are entangled states, one can find local operations L 2 
and -L3 acting on qubit 2 and 3 respectively such that 
|0) = Ii (g) L 2 \Y) and = I t <g> L 3 \Y) where \Y) is 
the sing let state (|01) - |10>)/v / 2- The ground space of 
IT12 + II13 is therefore 



L 3 \Y)(Y\ 13 I®Ll) 



K{I®L 2 \Y){Y\ 12 I 

= (it)-i/c(|y)(F| 12 + |F}(y| 13 ) 

= (L^- 1 fC(\Y){Y\ 12 + \Y)(Y\ 23 ) 
= /c(n 12 + L 2 ® L 3 \Y){Y\ 23 L\ ® 4), 

where the first equation uses Eq. (2), the second one 
is obtained by a direct calculation establishing that 
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/C(|F)(l r |i2 + |y)(i A |i3) is the symmetric subspace of the 
three qubits, and the last step employs Eq. (2) again. 
This validates the Type -I sliding operation. 

Repeated applications of the two types of sliding oper- 
ations can modify an arbitrary graph (a homogeneous 
Hamiltonian) with solid and dashed edges to the so- 
called simplified interaction graph (simplified homoge- 
neous Hamiltonian). The simplified graph has a back- 
bone of only dashed edges and several solid-edge tails 
attached to the backbone. An example of such a graph 
is shown in Fig. 3c. This simplification can be done in 
two steps by first changing each connected component of 
solid edges into a tail, and then sliding all dashed edges 
connected to a tail to one end of the tail. During the pro- 
cess of the sliding operations, it may happen that there 
is more than one edge between two vertices. If these mul- 
tiple edges represent different constraints, one will essen- 
tially have a high rank constraint and can deal with it as 
before in the case study of rank. 

Simplified homogeneous case. — Since sliding opera- 
tions do not change the ground space, we only need to 
work with simplified interaction graphs. The idea is to 
build the entire ground space by extending the ground 
space for the dashed backbone. Let us first consider 
the case where there is only one tail in the simplified 
interaction graph. More specifically, let S be the ground 
space of the dashed constraints in the backbone J, and T 
be the symmetric subspace confined by the tail of qubit 
set K, where J n K has exactly one qubit a, through 
which the tail is attached to the backbone. We prove 
that R = S (g> Hk\{o,} nT® %j\{ a } is again a product 
span. Write S as the direct sum 

(S ®H a )® (® S J a; „j , 

where |aj) Q 's are different dashed constraints on vertex 
a and d is number of such |<x,) a 's. For the basis of So, 
all the constraints in the backbone are already satisfied, 
and therefore, the qubit a can be any state. We say that 
qubit a is free in this case. For the basis of Sj, qubit a 
has to be \otj-) in order to satisfy all the constraints in the 
backbone. In this case, the state can only be extended 
to the tail by copying. In summary, the intersection R 
contains the space 

(S ®T)® (®Sj®\af)®\A . (4) 

We will need to show that this is actually everything in 
R. 

We first claim that the product basis for Sj's all to- 
gether form a linearly independent set. By orthogonal- 
ity (up to local operations). Sj and S& are orthogonal 
if \oij) and \ak) are not. On the other hand, if \aj) 



and \ak) are orthogonal, the basis for Sj and Sk are 
linearly independent. Otherwise, we will find a state 
in both Sj and Sk, meaning that \ip) should be in 
So, a contradiction. Now, for any state in R, we 
can write it as |\&) = J^j where |\£j-)'s are lin- 

early independent product states spanning S. Let \iffj) 
be the state on J\ {a} when the state on J is \^j). One 
can also collect terms according to the state on J\ {a}, 
that is, |*) = Ek\^k)Ei\^)t,i\^k,i). As shown pre- 
viously, |\I/fc)'s are linearly independent, and we know 
J2i \^)t.i\®k,i) is in T for each k. That is, the state \ip) 
is indeed in the space of Eq. (4) . As the symmetric sub- 
space can always be spanned by product states, we have 
finished the proof for the case of one tail. For multiple 
tails, the proof is essentially the same by an induction on 
the number of tails. 

Application to the counting of degeneracy. — The re- 
sults above actually allow us to prove that counting the 
ground state degeneracy of a 2BFF Hamiltonian is in $P. 
The class #P contains functions / if there is a polynomial 
time algorithm A such that 

f(x) = \{y,A(x,y) accepts.}), 

where y is usually called a proof to the verifier A. 

As indicated by the ground space structure in Fig. 1, 
the isometries will not change the dimension and we only 
need to consider the simplified homogeneous case where 
one can actually replace the solid edges of the tails to be 
dashed edges forming alternating chains. As long as we 
choose the constraint of the tail on the vertex connecting 
to the backbone to be different from all other constraints 
\aij) of that vertex, the dimension of the solution space 
remains unchanged. To understand this, we need to re- 
view the extension of the product span with intersection 
of symmetric subspaces. If the vertex in the intersection 
is free, we will have the whole symmetric subspace on 
the tail which is of dimension k + 1 where k the number 
of qubits in the tail. This coincides with the dimension 
of the alternating chain. If the vertex in the intersection 
is not free, we will have a unique extension in the tail, 
which again coincides with the case of alternating chain. 

It therefore suffices to count the dimension of any 
dashed graph. To show that it is in -jfP, one can choose 
the proof to the verifier to be the non-deterministic 0, 1 
choices in the case of (1) all-the-same-constraint vertex 
and (2) alternating loop. 
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